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We present a model which describes coherent and incoherent processes in continuous-variable 
atom-light interfaces. We assume Gaussian states for light and atoms and formulate the system 
dynamics in terms of first and second moments of the angular momentum operators. Spatial and 
temporal inhomogeneities in light and atom variables are incorporated by partitioning the system 
into small homogeneous segments. Furthermore, other experimental imperfections as for instance 
limited detector time-resolution and atomic motion are simulated. The model is capable of describing 
many experimental situations ranging from room temperature vapor cells to sub-mK atomic clouds. 
To illustrate the method, we calculate the effect of detector time-resolution, spatial inhomogeneities 
and atomic motion on the spin squeezing dynamics of rubidium 87 on the D2 transition. 
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I. INTRODUCTION 

In the last decade, many pioneering experiments 
have demonstrated quantum information processing with 
continuous variables in atomic and photonic systems. 
Atomic spin squeezing can overcome the standard quan- 
tum limit in magnetometry pQ. Macroscopic numbers 
of atoms can serve as a memory of a quantum state of 
light [2 . A variety of atomic systems have been pro- 
posed or demonstrated, including hot atoms in cells and 
cold atoms in magneto-optical traps and in optical dipole 
traps. The conditions vary greatly from system to sys- 
tem, notably in number of atoms, from ~ 10 12 in cells to 
~ 10 6 in dipole traps, in temperature, from ~ 300 K in 
cells to ~ 30 /xK in atom traps, and in the time-scale of 
the interactions, from milliseconds in cells to microsec- 
onds in cold atoms. At the same time, a variety of other 
physical effects, such as loss and decoherence of atoms, 
scattering and diffraction of light, and inhomogeneous 
coupling of the light and atomic variables are present to 
varying degrees in these systems. 

We demonstrate here that these many effects can be 
treated within a single framework. We work with Gaus- 
sian states, where coherent interactions and incoherent 
loss and decoherence processes have been studied. Previ- 
ous work on inhomogeneity has included mode matching 
[3] and introduction of weighted variables [J. In some 
cases decoherence effects due to inhomogeneous coupling 
have been identified j5j. Here we show how the model 
of Madsen et al. [6], when applied to the physical an- 
gular momentum and Stokes operators, can be naturally 
extended to include inhomogeneities as well as transport 
processes such as movement of atoms. 

In the first part we review some important definitions 
of angular momentum operators for atoms and light. We 
introduce the method of segmentation and give its math- 
ematical description in compact form. We give a general 



description of physical processes including coherent light- 
atom interaction, coherent and incoherent transport pro- 
cesses, and projective measurements. Incoherent trans- 
port is used to describe loss and decoherence. 

In the second part we apply these techniques to cal- 
culate the effect on spin squeezing of: imperfect detec- 
tor temporal resolution, spatial inhomogeneities in atoms 
and light, and atomic motion. 

II. PHYSICAL SYSTEM AND 
MATHEMATICAL MODEL 

A. Continuous variables for light and atoms 

Polarized light in the framework of continuous vari- 
ables can be described in terms of the Stokes operators 

w 

Here a = [&+ , a_] and a + , a_ are annihilation operators 
for circular plus and minus polarization, respectively and 
er x , (T y , and <t z are the Pauli matrices. The Stokes op- 
erators have the same commutation relations as angular 
momentum operators. In many situations of interest, one 
polarization component is strong. Here we consider lin- 
early polarized light: 

(S x ) = ^hN L = S x and ($„>,<&> =0 

where is the number of photons. A coherent polar- 
ization state can be expressed as an angular momentum 
- [3 or spin and atomic - j8] |9] coherent state. All of 
them have in common that the variances orthogonal to 
the main "spin" are 

var(S w ) = var(^) = hi 2 N h . 
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Along this direction we have either vai(S x ) = f^N^/A [7 
or var(S" x ) = jHE] . For AT L > 1, we can substitute the 
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operator S x by its expectation value. Hence, quantum 
polarization features are then solely contained in S y and 
S z . Geometrically, we are approximating a portion of the 
Poincare sphere as a plane, the geometry of the harmonic- 
oscillator phase space. Formally, this is referred to as 
the contraction from SU(2) to the Heisenberg-Weyl group 
|S]- Consequently, the commutator for S y and S z is not 
operator valued, as it would be in the SU(2) algebra. 
Instead of writing [Si, Sj] = ihsijkSk, we have [S y , S z ] = 
ihS x and [S x , S y / Z ] — 0. Apart from normalization, these 
are the commutation relations for the generators of the 
Lie algebra in the Heisenberg-Weyl group. 

For atoms, we similarly describe the collective spin of 
a collection of atoms with the angular momentum oper- 
ators 
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(2) 

defined in terms of bosonic operators b = [frj-, &jj . The 
states |t) , |4) are two degenerate atomic ground states. 
These could be the states of a spin- 1/2 atom, as in the 
proposal of Kuzmich et al. [10] , or more practically, two 
ground states of an alkali atom. Later, we will consider 
the case of F = 1, where |t) — \F = 1, mp = —1), |4) — 
\F = l,mp = +1) [TTJ. The operator J then describes a 
pseudo-spin, with angular momentum commutation rela- 
tions but without spin-like behavior under spatial rota- 
tions. 

We assume the atoms are polarized along a certain 
direction, so that one angular momentum component can 
be treated classically and the two orthogonal components 
carry the quantum properties. For x polarization, 
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and the variances are 



var( j y ) — var( j z ) — -h 2 N A 



B. Partitioning and covariance matrix 

A central goal of this work is to include spatial in- 
homogeneities in a description of the light-atom inter- 
action. In cell experiments, the atomic ensemble has a 
constant number density while a cold trapped sample can 
be highly inhomogeneous. In almost all experiments, the 
light distribution is inhomogeneous, e.g. from a Gaussian 
beam. 

We split the inhomogeneous ensembles of atoms 
and light into several segments. That is, we de- 
fine angular momentum variables for the atom seg- 
ments j( k,l \ and for the light segments S^ k ' 1 with 
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The segments orthogonal to the di- 
rection of light propagation are called transverse seg- 
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Figure 1: Schematic of partitioned atom-light interface. See 
text for details. 



ments (first index) and along the direction of light prop- 
agation, longitudinal segments (second index) (cf. Fig. 
[TJ. Here we do the segmentation in two dimensions only; 
the extension to the third is straightforward. In order to 
stay within the assumptions of the group contraction, we 
have to ensure that the particle number in each segment 



is itself large, i.e., N^ l > > 1 and A^ K '' ; > 1. The total 
angular momenta for atoms and light are 



(M') 



(k.i) 



and 



S = 



k.i 



k,V 



(3) 



We define a channel as the set of segments (light and 
atoms) which have the same transverse index k. 

We assume that the ensembles both of atoms and pho- 
tons can be described as Gaussian states in the harmonic 
oscillator phase space and that all operations we apply 
will map them into Gaussian states. For a single segment 
of light we can use equation ([T]) and define a phase space 
vector as 



3(fc.O = 
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and similarly for atoms 



(k,i) 
v 

(k,i) 



b (fe ' z) . (5) 



In a common phase space for atoms and light we define an 
overall phase space vector in terms of angular momentum 
operators as 



(6) 



which is readily rewritten as the direct sum of phase space 
vectors of the sub-systems for atoms (A) and light (L) 



(7) 



Gaussian states are completely characterized by their 
first and second moments. First moments (v) represent 
a displacement in phase space. Second moments or vari- 
ances are given by 
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|vAH(vA v) T |) - (v) A (v) 



1 /„ 
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(8) 



which is the covariance matrix. For our purpose of ex- 
amining entanglement and squeezing properties, only the 
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second moments are of interest. We can write the covari- 
ance matrix of the joint atom-light system as 
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A C 
C T L 



(9) 



where C describes correlation between atoms (A) and 
light (L). 



III. UNIFIED DESCRIPTION OF PHYSICAL 
PROCESSES 



character of J we take only magnetic fields along the z 
axis into account. Such a field results in a rotation about 
the z axis in the Bloch sphere. To ensure the validity of 
the group contraction we also limit the rotations to small 
angles. The Hamiltonian for the segment (k, I) is 



fr(k,l) 

magn 



fJsgF 



j(k,i) B (k,i) 



(13) 



Where, /is is the Bohr magnet on and gp the Lande 
factor. This description includes homogeneous as well as 
inhomogeneous magnetic fields. 



The dynamics of v and 7 is calculated by difference 
equations which describe small but finite changes be- 
tween time-steps. This allows us to model coherent light- 
atom interactions, losses and decoherence, measurement 
and transport processes in a consistent way. 

The phase space vector and covariance matrix are up- 
dated in finite time steps r as 



v(t- 



F r (v(t)) 



(10) 



The time step is chosen to be the duration of a longi- 
tudinal light segment. It should be short enough that 
a longitudinal segment can be considered homogeneous, 
but still contain many photons. 

We note that in most experiments the atomic sam- 
ple is much shorter than the coherence time of the light 
pulse, which means that in each channel only one longi- 
tudinal light segment will overlap with the ensemble at 
most times. In addition, the effect of the several longi- 
tudinal atomic segments is, from the light's perspective, 
sequential: s^™' 1 ) interacts with j^"' 1 ) then with j^™' 2 -*, 
and so forth. In the time t = (m — 1) r to t < mr, the 
gO,™) m teracts with all atomic segments in the n-th 
channel. 



A. Coherent effects 

For effects described by a Hamiltonian H which is lin- 
ear in the elements of the phase-space vector v, the phase 
space vector evolves as (to lowest order in r), 
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The last equality, which expresses the change in v in 
terms of a matrix T r , is possible by the linearity of the 
Hamiltonian and the c-number-valued commutation re- 
lations. The covariance matrix evolves as 



1 {t + r) = T Tl {t)T T T 



Single species effects 



(12) 



A magnetic field acts solely on the atomic spin and 
leaves the light unchanged. Owing to the pseudo-spin 



2. Atom-Light Interaction 

For a homogeneous system of light and atoms off- 
resonant interaction gives rise to an effective Hamilto- 
nian. For the F = 1 pseudo-spin system, it has the form 

(2) 

H cS <x (a' ' + ^-)5 io + a (1) S z J z + a^(Sj x + S y J y ) 

(14) 

where a^°\ and are the scalar, vector and tensor 
components of the polarizability |12j . 

For brevity, we will write this interaction as H c s(S, J). 
As the light pulse propagates through the medium, the ef- 
fects of ifefj^™' 1 ),^"' 1 )), fl e ff(s( n ' 1 \j( n ' 2 )), and so forth 
are applied in sequence to the covariance matrix. Note 
that loss and decoherence may be applied between these 
coherent evolutions. 



B. Noise considerations 

In addition to Hamiltonian evolution, loss, transport, 
and decoherence of atoms and/or photons can be de- 
scribed. These processes introduce extra noise into the 
system. A fully general description of a noisy Gaussian 
process is the Gaussian completely-positive map (GP), 
which acts on the covariance matrix as 



7' = M7M 5 



N 



(15) 



where the real matrix M transforms the phase space vec- 
tor and the real symmetric matrix N describes added 
noise. These must obey [13] 



N + iS' - iMSM T > 



(16) 



where illy = and S , similarly defined, are com- 

mutation matrices before and after the transformation 
(note that the commutation relations, which include the 
"classical" components J x , S x can change due to loss and 
decoherence). This places a lower limit on the noise in- 
troduced. Specifically, 



N = liS' -iMSM T | 



(17) 



where I • I indicates the matrix absolute value, is the min- 



imal symmetric matrix to satisfy ( 16 1 
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Loss and Decoherence from photon scattering 



Inevitably, the coherent interaction of equation (14 1 



will be accompanied by spontaneous emission of pho- 
tons, producing also incoherent changes in the atomic 
state. We use equation (151 to calculate the effect of loss 



and decoherence of atoms and photons. Here "loss" of 
atoms refers to the decay of atoms into meta-stable states 
which do not interact with the light. Decay of atoms into 
the |t), ||) states can cause decoherence of the spin state. 
While loss is not present in the ideal spin- 1/2 system pro- 
posed by Kuzmich et al. [10], in alkali metal atoms both 
processes are observed. For light there is no decoherence 
process since spontaneously emitted photons scatter into 
all possible spatial modes and are counted as losses. 
The covariance matrix transforms as 

7 (< + r) = M Tl (t)M? + N T , (18) 

where the decay is described by 

M T = (1 - rjr)Ia © (1 - e)h ■ (19) 

Here r\ T and e are scattering probabilities for an atom and 
a photon, respectively. For rubidium 87 these are given in 
terms of experimental parameters in the appendix. Noise 
will have the form N T = N T j oss + N T 
get 



with (17 1 we 



,loss 



ri T (l-Vr)^N A I 2 ®e(l-s)^N L h (20) 



and 



(21) 



Here p is the fraction of the scattered atoms which re- 
turn to the system, assumed to be in a mixed state. This 
model has been used in the literature |14] and serves to 
illustrate the method. A different model would be neces- 
sary to describe some processes, e.g., optical pumping. I2 
is the identity matrix in two dimensions. O2 is the zero 
matrix and reflects the fact that we don't consider any 
decoherence for the light. For all simulations that fol- 



low in section IV we assume we have exclusively atomic 
decoherence and no loss, i.e., p = l. 



C. Projective Measurement 

The next class of operations we can apply are measure- 
ments of atomic or light variables. While a measurement 
will collapse the value of an observable in a way that is 
fundamentally random, the resulting variances change in 
a way that is completely predictable: The variance of 
the measured observable becomes zero, the variance of 
the conjugate observable becomes large or infinite. The 
variances of other observables may also be reduced if they 
are correlated with the measured observable. 



A measurement can be described by a projection ma- 
trix V . For example, to measure a polarization compo- 
nent of the (n,i) light segment, Sg n ' 1 ^ = cos OSy™'^ + 
siaOSz = (Pg 1 '^) 1, ■ v, the projector would be the 



outer product V = Pg 1 ' 1 ' A p^™'^. In practical situa- 
tions, measurement of a light variable also implies that a 
light segment has reached a detector and thus is removed 
from the problem, reducing the dimension of the vector 
v. Upon measurement, the covariance matrix becomes 
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I7 - j(V~fV) j 1 



\(n,i) 



(22) 



Where |...|/ n ^ removes the column and row correspond- 
ing to the measured, and no longer existing, light segment 
(n, i). (■■■)" indicates the Moore-Penrose pseudoinverse. 



Equation (22 1 is well known in mathematical statistics 



to compute the conditional covariance matrix of multi- 
variate normal distributions [15] ■ A more detailed intro- 
duction of Gaussian operations on Gaussian states can 
be found in [TBI IT?] . 

For the calculations in Part 2, we consider a large-area 
detector, i.e., one which does not distinguish between 
different channels (see Fig. [TJ. Therefore, we define the 
measured light variable to be 



k=l 



(k.l) 



(23) 



D. Combining effects 

When several effects are present at the same time-step, 
they are applied sequentially. The order of application 
can influence the results of the calculation if the time- 
step is not small. For example, when the light-atom in- 
teraction of Eq. (12 1 and noise of Eq. (18 1 are both 



considered, we can have 



or 



7(t + t) = M r T r 7(t)T^M r + N 7 



7(i + r) = T r [M Tl (t)M r + N T ] T T T 



(24) 



(25) 



depending on which effect is applied first. Physically, this 
ordering has no meaning, and in the limit of small time 
steps t, both (24| and (25 I give the same result. In the 



simulations that follow, we reduce r until the effect of 
the ordering is negligible. 

IV. RESULTS 

Now we give three examples how the model can be 
applied in the context of atomic spin squeezing. For all 
simulations we consider a cold ensemble of rubidium 87 
atoms in a dipole trap. The set of used parameters can 
be found in the appendix. It is well known that for large 
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detunings from resonance we can reduce the dipole inter- 
action Hamiltonian ( 14 1 to 



time. Mathematically, the measured variable is the sum 
of all n light segments 



H, 



!ff (g(W),p.O) = fy^QWJM) 



(26) 



The coupling constant <jr( fc '0 is proportional to the vector 
polarizability and defined in the appendix. Note, the 
Hamiltonian does not explicitly depend on r because the 
Stokes operators are proportional to the flux of photons 
times t. As initial states we assume a pulse of hori- 
zontally polarized light, i.e., S x = NjJl/2 and a coherent 
superposition of the Zeeman substates \F = l,mp = —1) 
and \F = l,m F = 1), i.e., J x = N A h/2. 

The effect of the Hamiltonian p6| on the light is a 
rotation of S x about the z-axis by an amount propor- 
tional to J z . At the same time, J z is not altered in this 
process. A projective measurement of S y provides infor- 
mation about, and thus reduces the uncertainty of, J z . If 
both input states are minimum uncertainty states, spin 
squeezing is obtained. To monitor the evolution of this 
process we evaluate 2Var ( J z ) j J x which is also known as 
the spin squeezing parameter [18]. For squeezed states it 
will become less than unity. The smaller the spin squeez- 
ing parameter the higher the degree of spin squeezing. 
There are other criteria, for instance by Wineland et al. 
[1] derived in the context of precision spectroscopy. Re- 
gardless which of the definitions is applied, we obtain the 
same qualitative results. 

To make the comparison between different experimen- 
tal situations clearer, we normalize the timescale. We can 
define a time when the rotation of the light polarization 
due to the atom-light interaction (26 1 exceeds the shot 
noise of the photons, i.e., when the signal-to-noise ratio 
becomes one. We want this time to be characteristic for 
the system as a whole. Therefore, we neither partition 
atoms nor light and get 
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h 4 G 2 N A <f> 



(27) 



Where G is the collective interaction strength, N\ the 
number of atoms and 4> the photon-flux. A detailed 
derivation is given in the appendix. 



A. Detector time-resolution 

As a first example, we study the influence of the de- 
tector time-resolution on the amount of spin squeezing 
and show the importance of correct modeling of pulsed 
experiments even for pulses much shorter than the detec- 
tor time resolution. 

We define an ideal detector as one capable of detecting 
individual light segments. The covariance matrix would 



be updated in accordance to (22 1 each time a light seg- 



ment hits the detector. In contrast, we say a detector has 
no time-resolution if it detects all segments at the same 
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2 = 1 



and we apply the transformation ( 22 I to the whole covari 



ance matrix. (For simplicity we assume only one atomic 
segment. Nonetheless, we keep the transverse index to 
avoid confusion.) The projector has the form 



■P = <0> 2 ffi -U r , 
n 



(29) 



Where U n is the unit matrix of rank n. 

In Fig. [2] we show the results for both having (a) 
perfect and (b) no temporal resolution. In the case of 
no temporal resolution, the achievable spin squeezing is 
reduced at longer timescales. This is understood if we 
compare the information carried by different longitudi- 
nal light segments. Early light segments interact with 
the initial atomic state and later ones with a noisier ver- 
sion of it. If the detector is lacking temporal resolution 
all this different information is mixed. 

Now we compare the results to calculations which ne- 
glect all dynamics during the pulse duration, e.g., in 0]. 
We call this type of model "zero-dimensional" because it 
treats the light-atom interaction as a point-like event in 
time. Therefore, we assume that the light pulse is not 
partitioned into longitudinal segments. Curve (c) and 
(d) in Fig . [2] show the results if the noise is added after 
(cf. Eq. |24| and before (cf. Eq. 25 1 the interaction, 
respectively. It becomes obvious that even in the case of 
a pulsed experiment it is important to model light as a 
stream of sufficiently short segments. 



B. Spatial inhomogeneities 

In many experiments inhomogeneities in light or 
atomic distributions are present. We give two examples 
for typical situations that can arise. As the simplest test 
model we assume an atomic ensemble which consists of 
two equally sized transverse segments. 

For all the following calculations we assume an ideal 
detector. Furthermore, the total number of atoms Na, 
the photon flux <&, the total interaction cross section A, 
and all other parameters are fixed and stated in the ap- 
pendix. The figure of merit is the variance of J z for the 
complete atomic ensemble 



var ( J z ) = var 



(30) 



To verify the validity of the segmentation model, we 
consider first a homogeneous atom distribution either as 
a single or as two segments. In both cases we observe the 
same results in the presence of atom-light interaction, loss 
and decoherence and transport processes, independent on 
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Figure 2: The spin squeezing achieved by using an ideal de- 
tector (a) and a detector with no temporal resolution (b) are 
shown. The parameters are given in appendix A. Two zero- 
dimensional calculations are given for comparison. In curve 
(c) the noise due to decoherence and loss was added after and 
in (d) before the interaction. 



the segmentation. They reproduce the solid line (a) in 
Fig. §. 

The first example reflects the situation we would find 
for inhomogeneous light fields interacting with homoge- 
neously distributed atoms. We model this with two chan- 
nels of equal interaction cross-section A/2. We assume 
light is only present in one of the channels. The result 
is plotted as the dotted curve in Fig. [3] The overall 
spin squeezing is reduced. To explain this, we can eval- 
uate both channels independently. One channel contains 
all the photons and the maximal obtainable amount of 
squeezing will be the same as for the homogeneous distri- 
bution (solid line in Fig. [2j. This reflects a very impor- 
tant property in atomic spin squeezing. The achievable 
amount of squeezing does not depend on the intensity of 
light (supposed it is not zero) but rather on the optical 
depth of the atomic ensemble. For the second channel, 
without light, we expect no change in the atomic state. If 
we combine these two results, we get exactly the dotted 
curve shown in Fig. |3j 

The second example is the inverse situation. The light 
beam has a larger cross-section than the atomic ensem- 
ble. We model this case by assuming all Na atoms only in 
one of the channels and light homogeneously distributed 
over both. The result, plotted as the solid line in Fig. 
[3] seems surprising. We see the exact same dynamics 
as for the homogeneous case. In this situation two ef- 
fects are compensating each other. The optical depth for 
the atoms is twice as large as in the previous examples 
and leads to larger spin squeezing. On the other hand, 
the light which does not interact with the atoms is also 
detected, and contributes noise but no additional infor- 



mation about the atoms. 

The two examples give some intuition about the in- 
fluence of inhomogeneities. It is now straightforward to 
apply it to more interesting and complicated experimen- 
tal cases. 
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Figure 3: Two exemplary cases are compared to the homoge- 
neous situation were light and atoms are evenly distributed 
over space. For more information see text. 



C. Atomic motion 

As a last application, we ask what happens when atoms 
can change places and go from one segment to another. 
This is a relevant question comparing different experi- 
mental situations. Atoms in vapor cells, for instance, 
have approximately room temperature. This corresponds 
to a root-mean-square (rms) velocity of hundreds of me- 
ters per second as opposed to a few tens of millimeters 
per second for dipole-trapped atoms. Typical timescales, 
to, for the light-atom interaction are milliseconds and 
microseconds, respectively. Hence, atoms in vapor cells 
have moved around half a meter (effectively) whereas the 
atoms in the trap moved less than hundred nanometers. 
This suggests that for trapped atoms, any inhomogeneity 
in the light beam will be mapped onto them [23]. To find 
a more quantitative description, we introduce a mixing 
probability m T per time step r and per atom. It is de- 
fined as the probability an atom would escape from one 
segment to another in one time step. 

As in the previous section, we use the simple test model 
of two channels, where all the light is concentrated in one 
part and the atoms are homogeneously distributed over 
both. The real symmetric matrix describing the mixing 
is 



M = 



1 — m T 



1 — m T 



<h, 



(31) 
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and the introduced noise is given by 



CONCLUSION 



N = m T (1 - m T ) -^Na 



1 -1 
-1 1 



(32) 



The noise matrix N reflects two things. First, the in- 
dividual variance in a single segment will increase. This 
is not surprising, because already correlated atoms leave 
the segment and uncorrelated atoms enter. Second, the 
variances of the total spin components are not altered by 
mixing. This is expected, since the choice of the parti- 
tioning is arbitrary and can have therefore no influence. 
An alternate way to derive N and M is sketched in |2D 
As a concrete example we assume we have an ideal gas 
of atoms and derive the mixing probability from kinetic 
gas theory. The number of collisions per area and time 
in an ideal gas is known to be R = Nv lms / (V >/6n). 
Where, N is the number of atoms in volume V, and their 
rms velocity is v rms . From this we can calculate the rate 
at which individual atoms cross a surface of area A', r = 
RA' /N. Furthermore, we assume that the atoms occupy 
a box of volume V = A' \A~, where A is the interaction 
cross section of a segment. The rate can therefore be 
written as 



V6ttA 



(33) 



In the limit of small r we can define a mixing probability 
as m T = it. 



(34) 
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Where w rms = y/iksT ' jra is the rms velocity of the 
atoms. This model is not an exact treatment of the differ- 
ent physical situations we find in vapor cells and atomic 
traps. Nevertheless, it suggests how atomic motion influ- 
ences the formation of spin squeezing. 

In Fig. [4] we plot the squeezing factor for different 
mixing probabilities. For m r /r — * (blue curves) we 
have the same situation as in the dotted curve of Fig. [3] 
For increasing mixing probability (red curves) we see that 
the squeezing is improved and the full amount (compared 
to the homogeneous situation, i.e., solid curve of Fig. |2) 
is achieved again. In this limit, when the inverse mixing 
rate becomes the same order of magnitude as to, sufficient 
atomic movement is present that all the atoms get enough 
interaction to be uniformly squeezed. This suggest that 
approaches of matched variables, e.g., by Kuzmich et al. 
ID are more relevant for cold atoms than for hot vapor 
cell experiments. 

If we instead focus our attention only on the segment 
of atoms which is illuminated we see in plot b) of Fig. 
[4] that the spin squeezing for this segment reduces when 
the mixing probability increases. One can interpret this 
as a decoherence mechanism for the smaller segment |19] . 



We have presented a model to compute the dynamics 
of interacting light and atomic ensembles with Gaussian 
states. The model is based on covariance matrices for 
the quantum components of collective angular momen- 
tum operators and employs segmentation of the light and 
atom systems. The model is similar to that of Madsen 
et al. [6], but extends the segmentation to light and to 
the transverse directions. Also, we use angular momen- 
tum operators, rather than derived canonical operators, 
which give intuitive results and simplify partitioning. We 
show how to include many effects which arise in real ex- 
perimental situations, including spatial and temporal in- 
homogeneities, atom motion, loss, and noise introduced 
by photon scattering. 

Employing this model, we have made the following ob- 
servations: The dynamics of spin squeezing requires time- 
dependent modeling, even when the atoms interact with 
optical pulses which are shorter than the detection sys- 
tem can resolve. At the same time, the detector time 
resolution has only a minor effect on the degree of spin 
squeezing under realistic conditions. The effect of spatial 
inhomogeneities in light and atoms have non-equivalent 
effects on spin squeezing: Concentration of the light into 
a sub-region of the atoms produces equal squeezing of 
the sub-region, but less squeezing of the entire ensem- 
ble, while concentration of the atoms into one sub-region 
of the light gives equal squeezing of the spin ensemble. 
Finally, we observe that atomic motion between an illu- 
minated region and a non-illuminated one tends to de- 
grade squeezing of the illuminated region while increasing 
squeezing of the entire atomic ensemble. This suggests 
that high-fidelity experiments should use probe pulses 
which are either much shorter, or much longer, than the 
time-scale of the atomic motion. 

The model can straight-forwardly be adapted to more 
complicated experimental situations, for example a cold 
thermal cloud in a focused laser beam. Also, application 
to multi-pass schemes as proposed by by Takeuchi et al. 
|2Q] or Sherson et al. [22] is possible. 
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Appendix 

The atomic ensemble we consider has 10 6 rubidium 
87 atoms at a temperature of 30 /iK. For the light we 
assume a flux $ = 10 14 s _1 of linearly polarized photons 
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Figure 4: Atomic motion is simulated with different mixing rates, a) shows the degree of spin squeezing for the whole ensemble, 
b) for the illuminated and c) for the un-illuminated segment. The discussion is given in the text. 



with a detuning A = 1 GHz from the F = 1 -»■ F' = 
transition of the D 2 line. The corresponding wavelength 
in vacuum is A = 780.241 nm. Both atoms and photons 



which is valid for detunings much larger than the natural 
line width. 

We derive the characteristic time to- If we apply the 



The coupling constant g in (26 1 is directly related to interaction (26 1 for a time r we find 



the vector part of the polarizability tensor [21]. For the 
F = 1 hyperfine ground state in the limit of detunings 
larger than the natural linewidth we find 



9 



(*,0 



1 TX 2 1 



-4J (A)-5Ji(A)+5(J 2 (A)) 



Where A^ k,l> ) is the interaction cross-section of the seg- 
ment (k,l) and r/27r = 6.065 MHz is the natural line 
width of the 5P 3 / 2 excited state. The functions <5pv(A) = 
(A + Ao._F') 1 include the finite hyperfine splittings in 
the excited state: A 0i F' ls the hyperfine level spacing be- 
tween F' = and F' = 1, 2. It would be straightforward 
to include inhomogeneity in the local light shift, due to 
the dipole trap, as A^ fc '^. However, we don't consider it 
for the simulations given here. 

The scattering probabilities for photons and atoms are 
given by r) T — Ni, tT cr(A)/A and e = N A cr{A)/A, as in 
|14j . Where, A is the interaction area and Nz, iT and Na 
the number of photons and atoms in a segment, respec- 
tively. The off-resonant atomic scattering cross-section 
is 



a (A) = £S ( 4<5o(A)2 + 551 (A)2 + 762{A) ' 



(36) 



Sy 



Sy + GflJ z S X 



J y Jy + GtiS Z J X , 



(37) 



(35) where G = ^^(-iS (A) - 56, (A) + 5<5 2 (A)) simi- 



lar to (35 I. The variances are readily calculated for co- 



herent input states (unprimed) 



var(V) - 



-JVi 



L,t 



G 2 h 2 -N A -N, 
h 2 



L.Uj 



G A H A —N A N LM 



(38) 



The same also holds for var(j'). When the second term 
in brackets is unity this describes a signal-to-noise ratio 
of one and with A^l.t = &t this occurs for r = to as 



given in (27). If we use the parameters given above to is 
0.55 /its. 
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